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Abstract 

Charge-density-wave (CDW) order has long been interpreted as arising from a Fermi-surface 
instability in the parent metallic phase. While phonon-electron coupling has been suggested to 
influence the formation of CDW order in quasi-two-dimensional (quasi-2D) systems, the presumed 
dominant importance of Fermi-surface nesting remains largely unquestioned in quasi-ID systems. 
Here we show that phonon-electron coupling is also important for the CDW formation in a model 
quasi-ID system ZrTe3. Our joint experimental and computational study reveals that particular 
lattice vibrational patterns possess exceedingly strong coupling to the conduction electrons, and 
are directly linked to the lattice distortions associated with the CDW order. The dependence 
of the coupling matrix elements on electron momentum further dictates the opening of (partial) 
electronic gaps in the CDW phase. Since lattice distortions and electronic gaps are the defining 
signatures of CDW order, our result demonstrates that the conventional wisdom based on Fermi- 
surface geometry needs to be substantially supplemented by phonon-electron coupling even in the 
simplest quasi-ID case. As prerequisites for the CDW formation, the highly anisotropic electronic 
structure and strong phonon-electron coupling in ZrTe3 give rise to a distinct Raman scattering 
effect, namely, measured phonon linewidths depend on the direction of momentum transfer in the 
scattering process. 
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I. INTRODUCTION 


The mechanism for charge-density-wave (CDW) order in metals is a long-standing prob¬ 
lem in condensed matter physics. The recent discovery of ubiquitous charge-ordering phe¬ 
nomena in cuprate high-temperature superconductors [1-7] has aroused new interest in this 
problem, since a thorough and generic understanding of the formation of CDW order, even 
in simple metals, may shed light on several important issues relevant to the enigma of 
high-temperature superconductivity. In particular, the so-called pseudogap phenomena in 
underdoped cuprates [8, 9], which precede superconductivity upon cooling, resemble the 
opening of partial electronic energy gaps on certain parts of Fermi surfaces in CDW metals 
[10-12], While long-range static CDW order appears to compete with superconductivity 
[2, 5, 13], short-range CDW correlations are present over a substantial range of carrier con¬ 
centrations [2, 14] that support the appearance of superconductivity, which points to the 
intriguing possibility that the two phases might share a common microscopic origin. It is, 
therefore, important to identify the necessary and sufficient conditions for driving a CDW 
transition, and to understand how a material’s electronic structure can be affected upon 
approaching such transitions. 

CDW order is most commonly found in metals with highly anisotropic, or “low-dimensional”, 
electronic structure. Typical examples include quasi-one-dimensional (quasi-ID) K 2 Pt(CN) 4 
[15], blue bronzes [16], and transition-metal trichalcogenides [17, 18], as well as quasi-2D 
transition-metal dichalcogenides [19, 20] and rare-earth tritellurides [21]. In the classical 
explanation first introduced by Peierls [22] , CDW transitions are understood as arising from 
the presence of a pair of nearly parallel Fermi surfaces, known as Fermi-surface nesting 
(FSN), which causes the dielectric response of conduction electrons to diverge at the nesting 
wave vector q n . Even though a periodic modulation of the electron density can be stabilized 
only by coupling to the crystal lattice [23, 24], it is q n that determines the periodicity of 
the modulation, i.e., in the Peierls picture CDW order is primarily electronically driven. 

In most quasi-ID and some quasi-2D systems, a good agreement between the CDW wave 
vector q c dw and q n has indeed been found [12, 25, 26], in support of the Peierls picture. 

However, similar agreement between q c a w and q n is not supported by highly accurate 
measurements [10, 27, 28] and calculations [29] on the prototypical quasi-2D CDW compound 
NbSe 2 . First-principles calculations have cast further doubt on the generic validity of using 
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FSN as the sole criterion for potential Peierls instability [30]: on the one hand, FSN leads 
to a pronounced enhancement in the imaginary part of the electronic susceptibility near 
q n , but the formation of CDW order requires a large real part of the susceptibility which 
does not solely depend on FSN; on the other hand, the susceptibility enhancement due to 
FSN is usually weak in real materials, and it is excessively fragile against imperfect nesting, 
impurity scattering, and thermal-broadening effects. Analysis of FSN provides insights into 
the role of electronic instability in the formation of CDW order, but it is important to 
realize that the most predictive signature of potential CDW instability is the presence of 
soft phonon modes in the initiating phase [31]. In this line of thinking, it is only natural 
to recognize the essential role of phonon-electron coupling in the formation of CDW order, 
as suggested by several authors [30, 32], This idea has been explicitly tested in quasi-2D 
materials both theoretically [29, 33-36] and experimentally [36-38], and has proved quite 
successful in explaining the q cdw that deviates from q n . 
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FIG. 1. (a) Crystal structure of ZrTe 3 . (b) Temperature dependence of resistivity measured with 
currents flowing along the crystallographic a and b axes, normalized at 300 K. 

In quasi-ID materials, the importance of phonon-electron coupling relative to that of 
FSN has remained largely untested, as the good agreement between q n and q cdw hr ID 
cases might appear to render such tests unnecessary. One may note, however, that even 
in a simple quasi-ID system, the Peierls picture lacks such predicting power as to identify 
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lattice distortions that are associated with the CDW order. Moreover, CDW energy gaps 
are often found on only certain parts of the Fermi surface [12, 25], even when the remaining 
parts of the Fermi surface are equally-well nested, which cannot be explained by the Peierls 
picture. These inadequacies have motivated us to perform a detailed study of phonon- 
electron coupling in the prototypical quasi-ID CDW compound ZrTe 3 - ZrTe 3 possesses a 
monoclinic structure with a = 5.89 A, b = 3.93 A, c = 10.09 A, a = 7 = 90°, and f3 = 82.2° 
(Fig. 1(a)). CDW order is stabilized below T cdw = 63 K with q cdw = (0.07a*, 0,0.33c*) [39], 
which is in good agreement with q n that connects quasi-ID Fermi surfaces arising from the 
5 p orbitals of Te(2)/Te(3) atoms that form a chain-like structure along the a axis. 

ffere we combine variable-temperature Raman scattering measurements and first-principles 
calculations to investigate the role of phonon-electron coupling in the CDW order in ZrTe 3 . 
Despite the material’s simple crystal structure, the lattice vibrational properties of ZrTe 3 
have not been studied using first-principles calculations. A back-to-back comparison of 
our experimental and computational results yields excellent agreement, which allows us 
to reliably determine both phonon eigenvectors and phonon-electron coupling matrix ele¬ 
ments. We find that only phonons that involve longitudinal movement of Te atoms in the 
Te(2)-Te(3) chains exhibit large coupling to the conduction electrons. Importantly, we show 
that the electronic states that are most strongly coupled to the phonons are also those that 
become gapped out in the CDW phase. These results suggest that certain lattice vibrational 
patterns play a substantial role in the formation of the CDW order in ZrTe 3 , and that their 
interactions with the electrons dictate the opening of partial energy gaps on the Fermi 
surface. As a joint consequence of the highly anisotropic electronic structure and strong 
phonon-electron coupling, a distinct Raman scattering effect, namely, the measured phonon 
linewidths can depend on the direction of momentum transfer in the scattering process, is 
observed for the first time. 


II. METHODS 

A. Experimental methods 

Single crystals of ZrTe 3 were grown by a chemical vapor transport method using iodine as 
transport gas [40]. The crystals were characterized by x-ray powder diffraction and resistivity 
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measurements, which confirmed that our crystals were of single phase and high quality. The 
x-ray diffraction measurements were performed on a Rigaku MiniFlex diffractometer at 
room temperature. The resistivity measurement was performed with a standard four-probe 
method using a Quantum Design PPMS. 

Our Raman scattering measurements were performed in a confocal back-scattering ge¬ 
ometry on freshly cleaved crystal surfaces using the A = 514 nm line of an Ar laser for 
excitation. In order to obtain spectra with light polarizations along various crystallographic 
axes, crystals were cleaved both parallel and perpendicular to the ab plane. The Raman 
spectra were analyzed using a Horiba Jobin Yvon Lab RAM HR Evolution spectrometer, 
equipped with 1800 gr/mm gratings, a liquid-nitrogen-cooled CCD detector, and BragGrate 
notch Liters that allow for measurements down to low wave numbers. The temperature of 
the sample was controlled by a liquid-helium flow cryostat, with the sample kept under bet¬ 
ter than 5 x 10 - ' Torr vacuum at all times. Consistent Raman spectra have been obtained 
on several different samples, as well as with different excitation-photon wavelengths (633 
and 785 nm). 


B. Calculational methods 

Density-functional theory calculations were performed using Quantum Espresso [41], 
within the generalized-gradient approximation parameterized by Perdew, Burke and Ernz- 
erhof [42, 43], to investigate the electronic structure, BZ-center phonons and phonon- 
electron coupling in ZrTe 3 . Norm-conserving pseudopotentials, generated by the method of 
Goedecker, Hartwigsen, Hutter, and Teter [44] , were used to model the interactions between 
valence electrons and ionic cores of both Zr and Te atoms. The Kohn-Sham valence states 
were expanded in the plane wave basis set with a kinetic energy truncation at 150 Ry. The 
Raman-active phonons were calculated at the BZ center using a linear-response approach. 
The equilibrium crystal structure was determined by a conjugated-gradient relaxation, until 
the Hcllmann-Feynman force on each atom was less than 0.8 x 10 -4 eV/A and zero-stress 
tensor was obtained. A 12x18x8 k-grid centered at the T point was chosen, in combination 
with a Gaussian-type broadening of 0.0055 Ry. 

The phonon linewidth (defined as full width at half maximum, FWHM, in the Raman 
spectra), r y qu , is twice the imaginary part of the phonon self energy arising from phonon- 
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electron interactions: 
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where the summation is over a set of k-points forming a regular grid in the BZ, each 
with a pair of Kohn-Sham electronic states m and n whose crystal momenta differ by 
that of the phonon mode q. The phonon-electron coupling matrix elements, g^ n k + q m = 
(k.n\5V / + q m) / ^ 2cu qi ,, are obtained from Quantum Espresso [41], where u qi/ is the 

atomic displacements of mode v at a wave vector q. The Fermi-Dirac distribution, /k n , 
gives the occupation of the Kohn-Sham state with ek n - The summation in Eq.(l) can be 
quite singular, and we used the cubic-spline to interpolate the electronic spectra and the 
squared phonon-electron coupling matrix elements from a “coarse” k-grid of 44 x 68 x 24 to 
660 x 1020 x 360, where e kn were calculated. Finally, weighted average of 7 qt , from a small 
region of q near the BZ center was used to estimate the actual phonon line width measured 
in Raman experiments (see Sec. HIE). 


III. RESULTS AND DISCUSSIONS 

A. Sample characterization 

Figure 1(b) displays typical resistivity curves measured on our single crystals in the tem¬ 
perature range between 2 K and 300 K. A clear resistivity anomaly is observed around 
^cdw = 63 K when the electric current flows along the crystallographic a axis. In contrast, 
a similar resistivity anomaly is very weak or absent in the b axis. Filamentary supercon¬ 
ductivity is observed below T « 3 K. These results are consistent with previous reports 
[18, 45]. The small residual resistance at low temperatures demonstrates the high quality of 
our samples. 


B. Identification of Raman-active phonons 

The crystal structure of ZrTe 3 belongs to the space group P12i/ml, with all atoms 
occupying Wyckoff 2e positions. Group theory analysis shows that there are a total of twelve 
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Raman-active modes in ZrTe 3 , including eight A g modes and four B g modes, which involve 
atomic movements in the ac plane and along the b axis, respectively. Only linear photon 
polarizations are employed in our Raman scattering measurements, and the combination 
of incident- and scattered-photon polarizations is denoted by two letters that indicate the 
polarization directions with respect to the crystallographic axes. A g modes can be observed 
in aa, bb, and cc polarizations, whereas B g modes can be observed in cib and be polarizations. 




FIG. 2. Calculated vibrational patterns of eight A g (a-h) and four B g (i-1) phonon modes, from low 
to high energy, (m) Polarized Raman spectra obtained at 70 K, offset for clarity. Letters indicate 
the corresponding calculated vibrational patterns. 

Figure 2(m) displays Raman spectra obtained at T = 70 K > T c d w in the aa, cc, and be 
polarization geometries. The spectrum obtained in aa geometry is consistent with previous 
reports [46, 47] and reveals a total of six phonon peaks, but it fails to detect two remaining 
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Ag modes, presumably due to unfavorable Raman-scattering matrix elements. These two 
modes become clearly visible in the cc spectrum. The spectrum taken in be geometry reveals 
all four B g modes, which are found at different frequencies than those of the peaks in the aa 
and cc spectra. Cross-leakage signals related to the A g modes are present in the be spectrum, 
but are ignored in our identification of phonons. We note that some of the cross-leakage 
signals are clearly larger than expected if we only consider the imperfection of our polarizing 
optical set up (~ 1%), and they are hence likely related to the presence of disorder or defects 
in our crystals which may remove some of the symmetry selection rules. 

The successful observation of all twelve Raman-active phonons in Fig. 2(m) is crucial 
for the identification of their eigenvectors. Contrary to the assumption in previous reports 
[46, 47], the frequencies of the A g modes are not well-separated into a low- and a high- 
frequency group. Therefore the A g modes cannot be simply attributed to inter- and intra- 
ZrTe3-prism vibrations. We hence determined the phonon eigenvectors by first-principles 
calculations, the result of which is displayed in Figures 2(a)-2(l), where vibrational patterns 
are indicated by arrows with lengths proportional to the atomic displacements. Indeed, 
apart from a distinct difference between the A g and B g modes enforced by the structural 
symmetry, we find no clear regularity in the vibrational patterns. The calculated mode 
energies are displayed in Table I along with the measured values at T = 150 K. This tem¬ 
perature is chosen such that the phonons are least affected by anharmonic lattice interactions 
(stronger at higher temperatures) and the formation of CDW correlations (stronger at lower 
temperatures), which are not accounted for in our calculations. Despite a slight systematic 
underestimation in the calculated values by no more than 16 cm -1 , the overall agreement 
between the experimental and computational results is remarkable, rendering our determi¬ 
nation of the phonon eigenvectors highly reliable. A one-to-one correspondence between the 
calculated eigenvectors and the observed Raman peaks is shown in Fig. 2. 


C. Phonon linewidths 

An important observation in Fig. 2(m) is that the d and / phonon modes exhibit much 
larger linewidths than the rest. Moreover, the peak of mode d is distinctly asymmetric 
and can be better described by a Fano than a Lorentzian function. It therefore seems 
likely that these two phonons are strongly coupled to a continuum of excitations, such 
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TABLE I. Experimental and calculated phonon frequencies and linewidths at 150 K. All energy 
values are in units of cm -1 . The experimental determination of linewidth for modes e and g has 
large uncertainties because of their weak signal. 

as electronic excitations. To test this possibility we have performed variable-temperature 
Raman measurements (Fig. 3). By fitting the data to Lorentzian (mode /) and Fano (mode 
d) functions, we extract the evolution of peak parameters as a function of temperature 
(Fig. 4). The fitted parameters of mode h, which does not exhibit an unusual linewidth, are 
also shown in Fig. 4 for comparison. 

As temperature decreases, the linewidths (FWHM) of modes d and / are found to decrease 
substantially (Figs. 4(e) and 4(f)): the decrease starts already from above T cdw below a 
characteristic temperature T* m 140 K, and is particularly strong near T = T cdw . Meanwhile, 
as shown in Fig. 4(g) the asymmetry parameter q in the Fano function for mode d increases 
sharply in its magnitude as temperature is lowered below T cdw , which indicates that the mode 
rapidly recovers towards a Lorentzian lineshape as the CDW order develops. Importantly, 
we find that both modes d and / involve a longitudinal deformation of the Te(2)-Te(3) 
chains in their vibrational eigenvectors (Fig. 2), which will be discussed in more details 
later. For mode h, which does not exhibit a broad linewidth, the fitted FWHM evolves 
smoothly through T* and T cdw (Fig. 4(b)). Also the peak positions of modes d and / do not 
show clear anomalies near these temperatures (Fig. 4(h)). 

We attribute the higher characteristic temperature T* to the onset of CDW correlations 


9 

















FIG. 3. Raman spectra obtained with different photon polarizations at selected temperatures, 
offset for clarity. indicates phonon peaks that can be indexed by our calculational results 
shown in Fig. 2, indicates additional peaks that appear below T c dw, and “x” indicates signals 
related to cross leakage between different polarizations. 


in ZrTe 3 . A higher onset temperature of short-range correlations than the actual phase 
transition temperature can be generally expected in low-dimensional systems due to fluctu¬ 
ations. The decrease of phonon linewidths below T* thus suggests that the electronic states 
that are coupled to the phonons start to be removed from the Fermi surface by the incip¬ 
ient short-range CDW correlations. This is in qualitative agreement with ARPES results 
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FIG. 4. Temperature dependence of fitting parameters of selected phonons measured with aa 
polarizations. The peak indices ( d , /, and h ) are after those in Fig. 2. Vertical dashed line indicate 
Tcdw Shaded area indicates a higher characteristic temperature T* which can only be determined 
with a large uncertainty. In (a), (c), and (d), the peak areas are displayed after divided by the 
Bose factor. 

[12], in which T* is estimated to be above 200 K. In addition to the phonon linewidths, all 
the phonon intensities (integrated area under the peaks) show a clear departure from the 
high-temperature behavior below T* (Figs. 4(a), 4(c), and 4(d)), which suggests that the 
local structure of the crystal lattice starts to deform along with the development of CDW 
correlations. 


D. Electronic Raman signal 

To further examine the relationship between the opening of electronic energy gaps due 
to CDW correlations and the change in phonon-electron coupling, we first note that, es¬ 
pecially in the aa-polarized Raman spectra in Fig. 3(a), the background intensity exhibits 
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an anomalous decrease with decreasing temperature. This has also been pointed out in a 
recent report [47]. The highly accurate data in Fig. 3 allow us to fit all the phonon peaks and 
systematically remove them from the Raman spectra, a procedure that most clearly reveals 
the underlying electronic signals. Upon doing so, we find that there are further structures 
in the data, as shown in Fig. 5(a). Below ~ 100 cm -1 , the change of the electronic signal 
mainly occurs below T c d w , whereas the change between ~ 100 and 300 cm -1 occurs over a 
much wider temperature range, up to T*. This can be seen from integration of the electronic 
signal over different energy ranges (Fig. 5(b)). 



FIG. 5. (Color online) (a) aa-polarized Raman data obtained at selected temerpatures after sub¬ 
tracting all phonon peaks from the spectra, (b) Integrated intensities in (a) over different energy 
windows. Dashed line and shaded area indicate T c d w and T*, respectively, same as in Fig. 4. 

Being able to extract the electronic signal, we have attempted to correlate it with the 
change of phonon-electron coupling versus temperature. We first notice that the electronic 
signal (Fig. 5(b)) exhibits a similar T-dependence as the linewidths of phonons d and / 
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(Figs. 6(e) and 6(f)). Upon changing the energy-integration window (Fig. 5(b)), the degree 
of similarity varies, and the best agreement is found when we integrate the electronic signal 
up to 150 cm -1 . This “optimal” upper energy limit is in fact commensurate with the phonon 
energy scale. It implies that the phonons are damped by electronic excitations that are of 
lower energies than the phonon energies themselves. In addition to the linewidths, we find 
that the inverse of the Fano asymmetry parameter q for phonon d, which describes the 
peak’s departure from a Lorentzian lincshape due to coupling to an excitation continuum 
[48], is roughly proportional to the electronic signal’s amplitude over nearly a decade of 
change (Fig. 7). While this follows from the fact that the T dependence of q is not very 
different from that of the peak FWHM (Figs. 4(e) and 4(g)), it is reassuring to see that 1 jq 
extrapolates to zero in the limit of no electronic excitations - the peak would fully recover 
to a Lorentzian in the complete absence of phonon-electron coupling. 



FIG. 6. Temperature dependence of the aa-polarized electronic signal integrated between 20 and 
150 cm -1 , plotted together with that of the FWHMs of peaks d and / (Fig. 3(a)). Data are 
normalized at high temperatures. 


E. Momentum-resolved phonon-electron coupling matrix elements 

In order to reach a concrete and microscopic understanding of the experimental results 
presented in the previous Sections, we have computed the phonon-electron coupling matrix 
elements for phonons near the BZ center, as functions of the electronic band index and 
momentum k. Figure 8(b) displays our computed Fermi surfaces of ZrTe 3 , which consist 
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FIG. 7. The inverse of the Fano asymmetry parameter for phonon d versus the aa-polarized elec¬ 
tronic Raman signal averaged from 20 to 150 cm^ 1 , with temperature being the implicit parameter. 
Solid line is a linear fit to the data. 

of two main sectors as has been previously reported [12, 18]: A quasi-ID sector forms 
nearly parallel slabs running perpendicular to a*, which would indeed favor nesting with 
q n ~ q c dw = (0.07,0,0.33). The other sector is located primarily in the vicinity of the 
a* — c* plane. Four electronic bands contribute to the Fermi surfaces. They are marked with 
indices 1, 2, 3, and 4 in Fig. 8 (f), where their energy ranges are shown relative to the Fermi 
level. The quasi-ID sector of the Fermi surface located near the BZ boundary is derived 
from bands 3 and 4, which mainly arise from the 5 p orbitals of the Te(2)/Te(3) atoms that 
form a chain-like structure along the a axis (Fig. 1(a)). 

To reproduce the experimentally observed phonon linewidths, it is critical to realize that 
Raman scattering probes phonons in the small- but non-zero-q regime. In this regime, while 
the phonon-electron coupling matrix elements are well approximated by the q = 0 values, 
the summation of 7 in Eq. (1) differs dramatically from the q = 0 limit owing to intraband 
transitions. To account for this effect, we have weight-averaged 7 evaluated at a set of small 
q near q = 0 , with a maximum of \q\ = 6 x lCE 3 ^ 1 (ao = 0.5292 A is the Bohr radius) 
which corresponds to the back-scattering of 514 nm photons. An index of refraction of 4.6 
is used for calculating the momentum carried by 514 nm photons propagating inside ZrTe 3 
[49], and it is assumed that the matrix elements in this small region of q are the same as 
those at q — 0. For simplicity, the weight we used in the q-averaging is just 7 itself. This 
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FIG. 8. (a) BZ of ZrTe 3 . (b) Calculated Fermi surfaces, (f) is the corresponding energy locations of 
the bands across Fermi level (vertical line), marked by indices 1, 2, 3, and 4. (c)-(e) Isosurfaces of 
\dkn fen 1 2 f° r phonon modes d, /, and h in the BZ, with contributions from different bands relative 
to the chosen isovalues (vertical line) displayed in (g)-(i), respectively. The values of Ig^nk-rS 2 are 
magnified to similar amplitudes for easier visualization. The angle of view in (b)-(e) is similar to 
that in (a). 


way of weighing can be justified by the following argument. Both phonon-electron coupling 
and Raman process observe energy and momentum conservation, among other conservation 
laws. The computed phonon-electron matrix elements and hence y q naturally conforms 
to this constraint, and is therefore used as a first approximation to the likelihood of the 
particular electronic transitions involved in the Raman process. 

The calculated phonon lincwidths are listed in Table I along with experimental values. A 
semi-quantitative agreement is found between our experimental and computational results, 
and in particular we are able to obtain the three broadest peaks correctly in the computation, 
a result that does not depend on the detailed choice of the q region or the weighing function. 
This renders our determination of the phonon-electron coupling matrix elements highly 
reliable. 
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Indeed, we find that the computed phonon linewidths are dominated by intraband tran¬ 
sitions, especially within bands 3 and 4. An electronic wave vector k that makes large con¬ 
tributions to the linewidth should satisfy two conditions: (1) the matrix element | 2 is 
large, and (2) the electronic energy e*. n is sufficiently close to the Fermi level. To delineate 
the contributions to phonon linewidths from individual bands and different sectors of the 
Fermi surface, isosurfaces of intraband-transition contributions to the matrix elements for 
phonons d, /, and h are visualized in Figs. 8(c)-8(e), respectively. Figures 8(g)-8(i) display 
the corresponding value ranges of 1 9kn,kn\ 2 ar i s i n g from bands n = 1-4, relative to the iso¬ 
value for making the isosurfaces. The isosurfaces undergo continuous reduction when the 
isovalue is increased, and in particular for the phonons d (Fig. 8(c)) and / (Fig. 8(d)), they 
eventually disappear near the corners of the BZ. This shows, strikingly, that the phonons 
d and / are most significantly coupled to the electronic states near the boundary of the 
quasi-ID Fermi surfaces. It is not a coincidence that the CDW electronic gap is opened 
precisely in this region of the BZ [12], and we believe that it is this region of the BZ that 
makes the greatest contribution to the overall phonon-electron coupling, including to that 
of the acoustic phonon whose softening to zero frequency triggers the CDW transition (see 
Sec. Ill F) . The correspondence between the k regions shown in Figs. 8(c)-8(d) and the 
gap-opening regions reported in Ref. 12 suggests that the phonon-electron coupling matrix 
elements dictates the opening of electronic gaps on the Fermi surface. 

A remarkable commonality between phonons d and / is the involvement of longitudinal 
deformations of the Te(2)-Te(3) chains in their vibrational eigenvectors (Fig. 2). In contrast, 
phonons that do not involve similar atomic displacements all have narrow linewidths. This 
can be understood at an intuitive level, since it is the Te(2)-Te(3) chain that makes a major 
contribution to the electronic states near the Fermi level. The only exception is mode h, 
which also involves a deformation of the Te(2)-Te(3) chain but does not exhibit a very large 
linewidth in the Raman spectra. This is explained by the fact that the main contributions 
to ,9fc n kn 1 2 are from other regions of the BZ (Fig. 8(e)), and a considerable portion of those 
are not on the Fermi surface (Fig. 8(b)). It is hence no surprise that the linewidth of mode h 
does not show any pronounced anomaly near T cdw or T* (Fig. 4(b)); in contrast, a dramatic 
decrease of the linewidths and a recovery of a Lorentzian lineshape are found for mode d 
and / below T c d w (Fig. 3). 

A detailed comparison of the aa- and cc-polarization Raman spectra in Fig. 2(m) and 


16 


Fig. 3 suggests that the measured linewidth of phonon / is different between the two config¬ 
urations. In our backscattering experimental geometry, these two configurations correspond 
to the use of incident photons propagating along the c* and a* directions, respectively. This 
difference in linewidth is rather unexpected, and we are not aware of any previous report of 
similar effects. Our result suggests that the phonon linewidth sensitively depends on q very 
close to the BZ center. This dependence can actually be qualitatively reproduced in our 
computation: by choosing spherical q regions that are of the same size but centered around 
different positions relative to the T point for computing the weight-average of 7 , we find a 
50% larger linewidth for phonon / when the center resides along a* (Table I). A pictorial 
understanding of this result follows from the BZ regions visualized in Fig. 8 (d). The Fermi 
velocity vp in these regions are primarily along a*. According to Eq. (1), phonons with 
q vectors that satisfy q • vp ~ cu qi , will have the largest 7 due to intraband transitions. 
On the / phonon branch, we find that the cc-polarization measurement samples more of 
such phonons than the aa-polarization measurement, and hence observes a larger averaged 
linewidth. 

Despite the apparent similarity in the matrix elements of phonons d and / (Fig. 8 (c)- 
8 (d)), it turns out that the sampling for phonon d does not have a strong dependence on the 
center of the q region considered in our computation. While this result nicely corresponds to 
our experimental observations (Table I), it does depend on computational details including 
the precise band structure near the Fermi level, the shape and size of the sampling q region, 
and the weighing function, all of which are not known a priori. Therefore the agreement 
should be taken only at a semi-quantitative level. 

According to the above interpretation, the observation of different phonon linewidths with 
different incident-light directions requires the material system to have specific characteristics. 
First, the phonon linewidth must be large, and it must arise primarily from phonon-electron 
coupling that involves electronic intraband transitions (i.e., not from defects or anharmonic 
lattice interactions). Second, the main contributions to the coupling must come from a 
highly anisotropic part of the Fermi surface, such that the distribution of pertinent vp has a 
preferential direction. Failure to conform with the second requirement explains why phonon 
h does not exhibit a similar behavior as phonon / - the coupling matrix elements are large 
on the isotropic sectors of the Fermi surface (Figs. 8 (b) and 8 (e)). Finally, we note that 
the above two requirements are in perfect accordance with the characteristics of a “good” 
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CDW system, so it is not a coincidence that we can observe the effect in ZrTe 3 , which has 
arguably the simplest crystal structure (and hence phonon spectrum) among all quasi-ID 
CDW materials known to date. 


F. The CDW amplitude excitation 

So far we have only discussed q = 0 phonons which cannot be directly responsible for the 
CDW transition. In a recent x-ray scattering measurement of ZrTe 3 , a pronounced Kohn 
anomaly was found on the mostly transverse acoustic phonon branch with polarization 
along the a* direction [50]. The freezing of this mode to zero frequency was suggested 
to be the triggering factor for the formation of the CDW order. To further confirm this 
interpretation, we have searched for so-called CDW amplitude excitation in the A g Raman 
spectra [51]. Amplitude excitations can be generally expected in materials that undergoes 
incommensurate structural transitions including CDW and other forms of charge-ordering 
phase transitions [24, 51-53]. Such transitions can be generically understood as triggered by 
an anomalous phonon that freezes to zero frequency at the transition temperature. Deeply 
in the structurally distorted phase, because the structural instability becomes fully relaxed 
and/or the electronic states that exert damping on the anomalous phonon become mostly 
gapped, the frequency of the amplitude excitation is usually found to be close to that of the 
original (unrenormalized) phonon. In this regard, measurement of the amplitude-excitation 
frequency can be used to identify the key phonon. To the best of our knowledge, CDW 
amplitude excitation has not been reported for ZrTe 3 - 

Indeed, a close inspection of the low energy part of the Raman spectra in Figs. 3(a)-3(b) 
suggests the development of a new peak below 20 cm -1 at T = 10 K. This feature has been 
measured with very high counting statistics at low temperatures, as shown in Fig. 9(a), and 
above T = 45 K we can no longer resolve it above the diffused-scattering background. By 
fitting the peak to a Lorentzian profile, we have obtained its key parameters as functions 
of temperature (Fig. 9(b)-9(d)). In particular, its energy softens by over 10% upon heating 
from 10 K to 45 K. While such an amount of softening is considerably smaller than those 
observed in quasi-2D [19] and 3D materials [53], it is not uncommon among quasi-ID systems 
[54], in which the effect of fluctuations is expected to strongly reduce the actual transition 
temperature from its mean-held value. Therefore we attribute this feature to the CDW 
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amplitude excitation. A similar feature can be identified also in the spectra obtained with 
cc-polarizations (Fig. 3(b)), but it becomes very difficult to resolve above T = 10 K because 
of its weak intensity. Importantly, the amplitude-mode frequency (18 cm -1 ) at 10 K amounts 
to about 2.2 meV, not far from the unrenormalized energy (« 2.5 meV) of the Kohn anomaly 
at high temperatures [50]. This further confirms that, indeed, it is the acoustic phonon that 
involves the deformation of the Te(2)-Te(3) chains along the a axis that triggers the CDW 
transition at 63 K. 



12 14 16 18 20 22 0 10 20 30 40 50 

Raman shift ra (cm -1 ) Temperature (K) 

FIG. 9. (a) Raman spectra obtained with aa-polarizations at low temperatures, offset for clarity. 
(b)-(d) Temperature dependence of fitting parameters in (a). The peak areas in (b) are displayed 
after dividing by the Bose factor. 


IV. CONCLUDING REMARKS 

In summary, we have performed polarized Raman measurements on ZrTe 3 , the results 
of which are compared back-to-back with our first-principles calculations. The latter has 
allowed us to unequivocally determine the phonon eigenvectors as well as the phonon-electron 
coupling matrix elements, which have remained hitherto unknown for this model quasi- 
ID CDW compound. Our result demonstrates that particular lattice vibrational patterns, 
namely, the longitudinal deformations of the Te(2)-Te(3) chains, exhibit strong interactions 
with the conduction electrons. Such interactions might play an essential role in the formation 
of CDW order, even in our simple quasi-ID case where FSN has long been considered to be 
of dominant importance. 
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We emphasize that our calculations are for q « 0 phonons only. Our results about the 
phonon-mode- and clectron-k-dependence of the phonon-electron coupling matrix elements 
do not rely on the FSN geometry, but such dependences should be important to consider 
for phonons with q « q c dw as well, as they will cast a strong influence on the formation 
of CDW order. In particular, we believe that it is the matrix elements that select out 
the specific phonon that first freezes to zero frequency, and hence determine the lattice 
distortions that are associated with the CDW order. Moreover, it is highly likely that the 
opening of (partial) electronic gaps on the Fermi surface in the CDW phase is dictated 
by the k dependence of the matrix element for the selected phonons, especially when the 
quality of nesting does not vary substantially across a larger portion of the Fermi surface. 
Calculations for phonons with q q c dw warrant further investigations and appear necessary 
for a quantitative understanding of CDW order even in quasi-ID materials. 

As a joint consequence of the highly anisotropic electronic structure and strong phonon- 
electron coupling in ZrTe 3 , we have identified a distinct Raman scattering phenomenon, 
where the measured phonon linewidths can be pronouncedly dependent on the direction of 
the incident light, contrary to the usual understanding that Raman scattering is a q — 0 
probe. Similar effects can be expected to exist in other low-dimensional CDW materials, 
where the pertinent requirements are most likely to be conformed simultaneously. 
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